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NYO-1480-12 

Colllslonless  Shocks  and  Solitary  Waves 
CATHLEEN  S.  MORAWETZ 

Recently  there  has  been  a  revival  of  interest  in  the 
collision-free  magnet o-hydrodynamlc  shock.   This  interest  centers 
around  the  apparent  detached  shock  which  buffers  the  earth  and 
its  magnetic  field  from  the  solar  wind.   Observations  by  various 
satellites  indicate  that  a  phenomenon  occurs  as  in  Figure  1. 
Around  the  earth  there  is  a  large  vacuum  field.   A  stream  of 
plasma,  the  solar  wind,  has  to  pass  around  this  vacuum  field 
which  acts  roughly  like  a  body  in  a  supersonic  stream.   The  solar 
wind  is  slowed  down  by  crossing  a  shock  and  then  slips  around  the 
vacuum  field. 

The  Lundquist  equations  and  the  theory  of  the  detached  shock 

1  2 
corroborate  this  picture  qualitatively.  '    But  certain  details 

are  contradictory.   Namely  the  plasma  layer  between  the  shock 

and  the  vacuum  field  or  magnetosphere  appears  to  have  a  turbulent 

or  at  any  rate  highly  irregular  magnetic  field  and  plasma  velocity. 

Secondly  we  know  theoretically  that  there  are  virtually  no 

collisions,  hence  no  viscosity,  and  thus  the  classical  theory  to 


1.  John  R.  Spreiter  and  William  Pritchard  Jones,  Journal  of 
Geophysical  Research  68,  12  (1963)- 

2.  C.P.  Sonett,  Bull.  National  Academy  of  Sci.,  84  (l964). 
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justify  a  sharp  shock  Is  gone.   Finally  ahead  of  the  shock  edge 
the  orientation  of  the  field  has  an  oscillation. 

The  early  Interest  In  collision-free  shocks  was  mainly 
connected  to  the  problem  of  heating  rarefied  plasmas.    In 
ordinary  gases,  a  natural  method  for  friction  Is  to  compress  a 
gas  quickly  and  non-adlabatlcally  by  driving  a  sharp  shock 
through  It.   Can  the  same  effect  be  achieved  In  a  plasma  without 
standard  dlsslpatlve  forces? 

To  look  at  these  problems  theoretically  we  want  to  consider 
an  Idealized  mathematical  model.   The  plasma  Is  composed  of  two 
kinds  of  particles,  electrons  and  Ions.   We  suppose  that  there 
are  no  collisions  and  that  the  distribution  functions   f,(u,  x,  t) 
for  electrons  and  Ions  satisfy  the  Vlasov  equations.   The  electric 
and  magnetic  fields   (E,Bj  are  governed  by  Maxwell's  laws  where 
the  current   J  and  charge  density  q   are  found  as  Integrals  of 
the  distribution  functions   f   : 


C.S.  Gardner,  H.  Goertzel,  H.  Grad,  C.S.  Morawetz,  M.H.  Rose, 
and  H.  Rubin,  "Hydromagnetlc  Shock  Waves  In  High-temperature 
Plasma",  In  Proc.  of  the  Second  United  Nations  International 
Conference  on  the  Peaceful  Uses  of  Atomic  Energy,  Vol.  51, 
United  Nations,  Geneva  (1958),  230-237- 
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(1) 


(2)       V'B  =  0  curl  B  =  i^J 

kV'E  =  q  curl  E  +  ^  B  =  0 


(3)       J  =  e/(f^  -  f_)u|du| 
q  =  e/(f_^  -  f_)|du| 


Here   e   is  charge  per  particle,   m   Is  mass,   |j,  is 
permeability  and  K  is  small  in  our  units  of  length. 

We  consider  that  all  quantities  depend  on  only  one  space 
variable   x   and  time  t.   Furthermore  we  assume  charge 
neutrality,   (K  =  O),   instead  of  Poisson's  law: 

(^)    It  f^  +  u,  3-^  +  -  (E  +  uvB)-V  f^  =  0 
dt   +     1  dx     m  ^         '  ^  u  + 


(5)  It  f-  +  -1  ^  -  -f-  (E  +  uxB).V/_  =  0 

e  m 

^B  hB 

(6)  B^  =  constant,     -^  =    -\iJ^    ,  -^^  =   iiJ^    , 


With  three  space  variables  we  would  have  to  neglect  second 
moments  involving  magnetic  fields  and  at  most  one  flow 
quantity.   For  high  pressure  plasmas  this  is  more  reasonable, 
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as,     bE  SB       aE 

(^)    ^  =  0  '  ^- ar  =  0  '  at^+ al^=  0 


Here   e  =  m_/m   . 

The  first  standard  step  is  to  take  moments.   We  denote  by  <(Q)  , 
the  integral   /Q,|du|.   The  equations  for  the  conservation  of  total 
mass,  total  momentum  and  total  energy  are  then: 


On  the  other  hand  from  the  Ion  equation  alone  we  shall  use 
only  conservation  of  mass  and  momentum  in  the   y,z   directions: 

(ii)|t<^>-^<V.)=o 

(12)  1^  /uf  N  +  1^^  4'l''^+)  "  ^  V^  ^  uxB)f  \  =0      (y,z   components) 

Next  we  integrate  these  equations  and  Maxwell's  equations 
with  respect  to   t,   average  over  a  large  time   T  and  assume  that 
there  are  no  gross  changes  with  respect  to  time,  i.e.. 


+T 

m    ^  F  dt  =  0.   Such  states  will  be  called  quasi-steady.   We 

J  _T 
are  left  then  with  steady  state  equations  for  the  time  average 

moments  and  the  time  average  fields.   With  some  computation  and 

Integrating  we  obtain  the  total  conservation  equations: 

(13)  Mass  <(u(f^+e^f_)y  =  pu  =  const. 

(14)  Mom.  /u  u(f  +e^f_)^  +  ^i 


1  '  +     -  V   m_|_ 


i(B^+Bp 


B  B   ) =  const, 


V   ^1^3 


(15)  Energy  /|u2 1  u^(f_^+e^f_)\  +  |-  (E^B^  -  E^B^  )  =  const, 


(16)  E  =  E  =  const.      B^  =  const. 

where  all  quantities  are  time  averages. 

From  the  Ion  equation  (4)  by  taking  moments  we  obtain 
similarly  ^(^^1^+^="   0   °^ 


(17)  <^-i^4.>=  constant 

(18)  |^<u,uf^y^<^(E  +  u  B)f^^ 

We  finally  average  the  last  equation  by  considering 


-  0. 


ex 


2 
dx    operating  on  it.   Neglecting  terms  of  order   e   and 


-X 

using  the   conservation   laws    (l6)   and    (l?)    one   finally   obtains 
E  +  uxB=0        in      y,z      directions. 

Here  the  "bar  denotes  the  average  over  both  space  and  time, 
i.e.: 

-  ^  fuf [duldxdt 


u 


/f [duldxdt 


We  also  take  space  averages  on  the  conservation  laws  and  we 
obtain  finally  the  transition  conditions  for  quasi-steady  states 
which  are  the  same  as  ( 1^-16)  with  space  and  time  averages  Instead 
of  time  averages.   Introducing  the  pressure  tensor  and  assuming 
that   /  f  dxdt   is  symmetric  about  the  mean  velocities-^  we  find  next 
the  transition  relations  are  the  deHoffman-Teller  relations  for  a 
gas  with   Y  =  -^  .   This  means  that  even  for  a  turbulent  flow  the 

standard  shock  and  contact  discontinuity  relations  hold   if  the 

3 
fine  structure  is  small  enough. 


5.  Note  we  have  assumed  no  dependence  on  y  and   z. 

6.  One  might  ask  where  is  the  corresponding  phenomenon,  in  fluid 
dynamics^  but  this  probably  does  not  occur  because  the  laminar 
flow  ahead  of  the  shock  has  a  higher  Reynolds  number  than 
the  flow  behind  It. 
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It  is  therefore  not  surprising  that  the  detached  shock  data 
computed  by  Sprelter  and  Jones  (see  footnote  l)  fit   quite  well. 
The  major  error  occurs  in  computing  the  temperature,  i.e.,  in 
one  of  the  higher  moments. 

The  next  problem  is  to  examine  the  fine  structure.   Only  the 
zero  temperature  case  has  been  completely  investigated  and  is  truly 
valid.   But  various  pressure  tensor  hypotheses  confirm  that 
qualitatively  in  most  cases  it  is  right.   In  this  case  the  Vlasov 
equations  reduce  to  six  ordinary  differential  equations,  the 
equations  of  motion  of  ion  and  electron.   We  first  ask  what  steady 
solutions  of  the  one-dimensional  equations  exist.   A  major 

qualitative  distinction  must  be  made  between  situations  where  there 

7-13 
is  and  where  there  is  not  a  magnetic  field     in  the  direction 

of  propagation;  case  i,   B-,  =  0  and  case  ii,   B.,  ^   0. 

7.  K.W.  Morton,  J.   Fluid  Mech.  l4,  part  3,  369-384  (1962). 

8.  A.  Banos,  and  A.R.  Vernon,  Nuovo  Cimento  XV,  269  (1960). 

9.  J.H.  Adlam  and  J.E.  Allen,  Phil.  Mag.  3,  ^^8  (1958). 

10.  K.W.  Morton,  Phys.  Fluids  J^,    I8OO  (l964);  "Finite  Amplitude 
Compression  Waves  in  a  Collision-Free  Plasma",  NYO-10434, 
MF-36,  Courant  Inst,  of  Math.  Sci. ,  New  York  Univ.  (l964). 

11.  P.G.  Saffman,  J.  Fluid  Mech.  1^,  552  (1961). 

12.  R.  Peyret,  C.R.  Acad,  Sc .  Paris,  t.258  ( I6  mars  1964) 
Groupe  2,  2973-2976. 

13.  A.D.  Pataraya,  Zh.  T.F,  32,  5  (1962)  (Soviet  Phys.,  Tech. 
Phys.  7,  5  (1962). 


Two  possibilities  occur.   Either  there  exists  a  solitary 
wave  and  a  family  of  periodic  solutions  or  there  exist  only  a 
constant  state.   For  a  fixed  angle  between  the  magnetic  field 
and  the  direction  of  propagation  the  solitary  wave  exists  if 
the  speed  of  propagation  is  small  enough. 

If  there  exists  a  solitary  wave  and  if  the  equations  of 
motion  are  modified  to  include  an  arbitrarily  small  amount  of 
friction  Morawetz   for  case  (i)  and  Morton  (see  reference  10 ) 
for  case  (ll)  have  shown  that  there  exists  a  steady  state  shock 
structure;  i.e.  there  is  a  solution  leading  from  one  constant 
state  to  another.   This  structure  has  the  following  properties: 

(i)   If   B  =  0,   there  is  a  front  propagating  with  the 
usual  shock  speed  followed  by  an  infinite  train  of  oscillations 

in  all  quantities  with  a  characteristic  wave  length  equal  to  the 

IS 
geometric  mean  of  the  two  phase  lengths.  -^     The  damping  of  the 

oscillations  tends  to  zero  as  friction  vanishes. 

(ll)   If   B  7^  0,   there  is  a  front  propagating  with  the 

fast  shock  speed  and  preceded  by  precursors  oscillations  in  the 

field  orientation  only.   Two  characteristic  wave  lengths  occur, 

the  ion  and  electron  phase  lengths.   In  the  limit  of  zero  mass 

ratio  only  the  ion  phase  length  occurs. 


\K.      C.S.  Morawetz,  "Magneto-hydrodynamlc  Shock  Structure  Using 
Friction",  NYO-8677,  Inst,  of  Math.  Sci.,  New  York  Univ. 
(1959). 

15.   The  distance  a  particle  travels  along  a  magnetic  line  during 
a  complete  change  phase. 
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If,    on  the  other  hand,  the  solitary  wave  does  not  exist  then 
neither  does  a  shock  structure  exist.   Morton  (see  also  reference 
16),  has  computed  the  solution  of  a  time-dependent  problem.   An 
Imposed  magnetic  field,  see  Figure  2,    drives  a  plasma  ahead  of  a 
vacuum  field.   If  the  imposed  magnetic  field  is  not  too  large  there 
is  a  continuous  solution  which  at  a  fixed  large  time  looks  like  the 
corresponding  friction  solution  described  above. 

If  the  imposed  field  Is  large  then  there  does  not  exist  a 
continuous  solution.   Morton  admits  an  inner  shock  connecting  two 
smooth  states  but  Its  structure  is  unknown.   Auer  et  al.    find  a 
turbulent  flow  for  the  case  B-,  =  0  but  it  is  not  clear  that  any 
actual  turbulence  will  have  this  particular  fine  structure  since 
there  would  surely  be  velocity  components  in  the  y,z  directions 
for  most  particles. 

Even  in  the  case  of  existence,  turbulence  may  be  caused  by 
Instability.   If  B  =  0,   the  charge  separation  field  E  is 

large  for  small  mass  ratio.   The  electrons  drift  rapidly  across 

17 
the  Ions  and,  as  was  pointed  out  by  Kellog,    the  whole  flow 

probably  exhibits  two-stream  instability. 

If  B  7^  0  we  might  also  expect  instability  in  the  front. 

The  precursor  orientation  oscillation  may  destroy  the  symmetries 

"I  R 
in  the  pressure  tensor  and  Llist's  conditions  for  stability   be 


16.  P.L.  Auer,  H.  Hurwltz,  Jr.,  and  R.W.  Kilb,  Phys .  Fluids  ^,    3 
(1962). 

17.  P.J.  Kellog,  Phys.  Fluids  7.  1555-1571  (1964). 

18.  R.  Liist,  "On  the  Stability  of  a  Homogeneous  Plasma  with  Non- 
Isotropic  Pressure",  in  U.S.  Atomic  Energy  Report  TID-7582, 
154-157  (1959). 
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violated. 

If  we  examine  the  data  (see  reference  2)  which  are  available 
for  the  IMP  satellite  we  see  that  generally  the  shock  will  have 
B  7^  0  and  a  study  of  the  field  shows  that  some  kind  of  precursor 
ahead  of  the  shock  edge  occurs.   In  fact,  the  state  ahead  is 
roughly  constant  but  for  an  oscillation  in  the  orientation.   The 
data  is  unfortunately  not  such  that  we  can  find  the  characteristic 
wave  length  or  even  tell  whether  it  exists  at  all. 

As  for  the  problem  of  heating  a  plasma,  we  can  only  say  that 
turbulent  shocks  should  probably  be  expectedo   The  desired 
temperatures  may  not  be  realized  because  much  energy  is  in  the 

turbulence.   But  nevertheless  the  detached  shock  due  to  a  magnetic 

19  /  ^ 

singularity    (for  the  earth,  a  dipolej  suggests  a  mechanism  for 

ionizing  and  heating  a  plasma  v/hich  would  then  be  confined  partly 

by  a  concave  stable  front  with  a  vacuum,  partly  by  a  shock  front, 

and  otherwise  by  the  unionized  gas  stream.ing  by. 


19.   J.  Hurley,  Phys.  Fluids  4,  109  and  854  (1961) 
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